Theory for superconductivity in a magnetic field: A local approximation approach 
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We present a microscopic theory for superconductivity in a magnetic field based on a local approx- 
imation approach. We derive an expression for free energy density _F as a function of temperature 
T and vector potential a, and two basic equations of the theory: the first is an implicit solution 
for energy gap parameter amplitude |Ak| as a function of wave vector k, temperature T and vector 
potential a; and the second is a London-like relation between electrical current density j and vector 
potential a, with an "effective superconducting electron density" Ua that is both T- and a-dependent. 
The two equations allow determination of spatial variations of a and jAkj in a superconductor for 
given temperature T, applied magnetic field Ha and sample geometry. The theory shows the exis- 
tence of a "partly-paired state," in which paired electrons (having |Ak| > 0) and de-paired electrons 
(having ]Ak| = 0) co-exist. Such a "partly-paired state" exists even at T = when Ha is above a 
threshold for a given sample, giving rise to a non-vanishing Knight shift at T = for Ha above the 
threshold. We expect the theory to be valid for highly-local superconductors for all temperatures 
and magnetic fields below the superconducting transition. In the low-field limit, the theory reduces 
to the local-limit result of BCS. As examples, we apply the theory to the case of a semi-infinite 
superconductor in an applied magnetic field Ha parallel to the surface of the superconductor and 
the case of an isolated vortex in an infinite superconductor, and determine, in each case, spatial vari- 
ations of quantities such as a and |Ak|. We also calculate magnetic field penetration depth A(r, Ha) 
and lower critical magnetic field Hci{T). The ratio Hc\{T) / HdT) (where He is the thermodynamic 
critical magnetic field) is found to be only weakly T-dependent for low temperatures and nearly 
T-independent for intermediate and high temperatures, and quantitatively not very different from 
that of the Ginzburg-Landau theory for Ginzburg-Landau parameter «: ^ 1. 

PACS numbers: 74.20.-z, 74.20.Fg 
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I. INTRODUCTION 

Currently accepted microscopic theories for supercon- 
ductivity in a magnetic field include the microscopic 
derivation of the Meissner effect in the BCS theorjii 
and the microscopic derivation of the phenomenological 
Ginzburg-Landau equations by Gorkov^ from the BCS 
theory. However, the BCS derivation of the Meissner ef- 
fect is based on a linear-response approach (i.e., an exter- 
nally applied magnetic field is treated as a weak pertur- 
bation), and thus, is valid only in the low-field limit. The 
Gorkov derivation of the Ginzburg-Landau equations is 
based on the assumption that energy gap function A(x) is 
a small quantity, and thus, is valid only for temperatures 
just below critical temperature T^- Various extensions of 
the Gorkov theory also rely heavily on the assumption 
of A(x) being small, or the assumption of magnetic field 
being weaki^ 

It is desirable to have a microscopic theory for su- 
perconductivity in a magnetic field that is valid under 
more general conditions. As an effort along this line, we 
present in this paper a theory based on a local approxima- 
tion approach. The theory allows microscopic description 
of the suppression of superconductivity by an externally 
applied magnetic field. The theory is expected to be valid 
for highly-local superconductors (the high-temperature 
copper-oxide superconductors arc examples of highly- 
local superconductors, for which Ginzburg-Landau pa- 
rameter K ^ 1) for all temperatures and magnetic fields 
below the superconducting transition (except, perhaps. 



in the high-field limit when the effect of spin param- 
agnetism, which is neglected in the present work, may 
become important). In the low-field limit, the theory 
reduces to the local-limit result of BCS.^ 

In Sec. m we present the details of the theory, and de- 
rive an expression for free energy density as a function 
of temperature T and vector potential a, and two basic 
equations of the theory: the first is an implicit solution 
for energy gap parameter amplitude |Ak| as a function 
of wave vector k, temperature T and vector potential a; 
and the second is a relation between electrical current 
density j and vector potential a. We also analyze solu- 
tions for |Ak(T, a)|. In Sec. IIIIl as examples, we apply 
the theory to the case of a semi-infinite superconductor 
in an applied magnetic field Hq parallel to the surface of 
the superconductor and the case of an isolated vortex in 
an infinite superconductor, and determine, in each case, 
spatial variations of quantities such as a and |Ak|. We 
also calculate magnetic field penetration depth A(T, Ha) 
and lower critical magnetic field Hci{T). A brief sum- 
mary is given in Sec. IIVI 



II. THEORY 

We consider a superconductor in an applied mag- 
netic field. Our starting point is the same BCS pairing 
Hamiltoniani 

H =Xl(^k-eF)4^Cka+^ Vkk'4^cli^;C_k'iCk'T> (1) 
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kk' 
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where ek is the normal state single-electron energy, 
the Fermi energy, y^k/ the pairing interaction matrix 
element, and c\^^ and Cko- are the Fermi operators of 
an electronic state of wave vector k and spin a in the 
normal state. Single-electron energy Ck and its corre- 
sponding single-electron wave function (/>k(x) satisfy the 
Schrodinger equation 



Since the non-local effect (or coherence effect in the 
superconducting state is not accounted for, we expect 
this approach to be valid only for highly-local supercon- 
ductors for which magnetic field penetration depth A is 
much larger than coherence length ^ (i.e., A 3> ^). 

We explain the details of this approach in the following 
subsections. 



if(x)0k(x) = ek0k(x) 
with single-electron Hamiltonian 



i/(x) 



1 
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-iW + -a(x) + C/(x), 



(2) 



(3) 



B. 



where a(x) is the vector potential, and U (x) a periodic 
scalar potential. 

In writing down Hamiltonian _ff , we have neglected, for 
simplicity, the effect of spin paramagnetism (which may 
be important for high fields). 

We note that, if ek and (?!)k(x) for a ^ are known, 
H can be diagonalized in essentially the same way as 
for a = (i.e, by making the BCS pairing approxima- 
tion, and then applying the Bogoliubov transformation*). 
However, since a(x) in the superconducting state is itself 
an unknown function, a simultaneous determination of 
Ek, '/'k(x) and a(x) does not seem possible. We therefore 
adopt a local approximation approach, which we describe 
next. 



A. Local Approximation Approach 

We note that, in the superconducting state, vector po- 
tential a(x) varies spatially with the length A, the mag- 
netic field penetration depth, which is '-^ lO^A. In con- 
trast, single-electron wave function (/'k(x) oscillates spa- 
tially with a much shorter length kp^^ lA (here \s.f is 
a Fermi wave vector) so that a(x) can be considered 
locally constant over many wavelengths of (/)k(x). Our 
approach is based on this observation, and can be out- 
lined as consisting the following three steps. 

Step 1: We focus on a small region around a local 
point, say x'. Dimension D of this small region satisfies 
d D X, where d ~ lA is a crystal lattice constant 
of the superconductor. In this small region, we assume 
a(x) = a(x') is a constant, and solve the Schrodinger 
equation [Eq. ([2])] to obtain ek and (/)k(x). 

Step 2: Based on the obtained ek and 0k (x) in the 
vicinity of x', we diagonalize Hamiltonian iJ, which now 
becomes a local quantity, because of its dependence upon 
a(x') through ek and 0k (x). 

Step 3: Once Hamiltonian H is diagonalized, we move 
on to derive an expression for local free energy density 
F and two basic equations of the theory. The equations 
allow determination of the spatial variations of vector 
potential a and other quantities for given temperature, 
applied magnetic field and sample geometry. 



Determination of ek and (;/)k(x) 



As outlined above, we first focus on a small region 
around a local point, say x'. Dimension D of this 
small region satisfies d ^ D ^ X, where d is a crys- 
tal lattice constant of the superconductor, and A is the 
magnetic field penetration depth. Since vector poten- 
tial a(x) varies spatially with the length A, which is 
^ 10'^ A, whereas 0k (x) oscillates spatially with the length 
d, which is ^ lA <C A, we can assume a(x) = a(x') ~ 
constant in this small region, and solve Eq. ([2]). 

For a constant a(x) = a(x'), it is not difficult to solve 
Eq. (21). For free electrons, for which scalar potential 
U{x.) — constant, solutions of Eq. ([2]) are easily obtained. 
Namely, for a = 0, we have 



.(0) 



and 



e^ ' — h^k'^ /2m 



and for a 7^ 0, we have 



Ck — (jik. H — a^ /2m 



and 



0k (x) = 



Here we have used ej^^"* and 0j^"'' to denote solutions of 
Eq. dal) for a = 0. 

In this paper, we will not consider the case of a general 
periodic scalar potential C/(x) (i.e., we will not consider in 
this paper how the details of an electronic energy band 
structure may affect properties of the superconducting 
state). Instead, for simplicity in presenting the theory, 
we will use the solutions for free electrons in the following. 



(4) 
(5) 

(6) 
(7) 



C. Diagonalization of H 

Having obtained ek and 0k (x), we can move on to diag- 
onalize Hamiltonian H of Eq. ([T]), which now becomes a 
local quantity, because ek and 0k (x) are obtained locally 
at x' for a — a(x'). 

An important step in the diagonalization of Hamilto- 
nian H is assuming that the pairing^ of electrons of op- 
posite momenta and spins holds even for a ^ 0. I.e., for 
a pair of (kt) and (— k|) electrons, we assume 



-kiCkt) 7^ 



(8) 
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when the electrons are superconducting (here the angle 
brackets (• • •) denote a thermal average). 

The diagonalization of Hamiltonian H is the same as 
in the case of a = except that we have e_k ^ eic for 
a ^ 0. The results of the diagonalization are as follows. 

The energy gap parameter is defined aa^ 



Ak = - ^ Vk,k' (c-k'xck't) 

k' 



(9) 



The diagonalized Hamiltonian is 
H = Y.(U^ + i?k7,^T^kT + i?-k7lki7-ki ) , (10) 



where 



C/k = ^^:^-4^^ + ^(l-/k-/-k); (11) 



(12) 



is the single-electron energy in the normal state, mea- 
sured relative to the Fermi energy e^; 



£^k - 



(s) ^ Ck - C-k 
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(13) 



the quasi-particle excitation energy in the superconduct- 
ing state; 



the symmetric part of _Bk; 



/k= (e^-/'=-^ + l 



(14) 



(15) 



the Fermi function; and 7][.^ and 7k(T are the Fermi oper- 
ators for quasi-particles in the superconducting state. 

The operators 7^^ and 7kcr are related to cj^.^ and Cko- 
via the Bogoliubov transformation^ 



CkT 



"k 



-f'k "k y V "^-ki 



7kT 



(16) 



where the coefficients Mk and Uk satisfy the following re- 
lations: 



(17) 




2E, 



and 



AkWk^k 



2E, 



is)- 



(18) 



(19) 



After the diagonalization of Hamiltonian H, Eq. Q 
can be expressed as 



1 - /k' - /-K 



-Ak 



(20) 



With regard to the above-described diagonalization of 
Hamiltonian H, the following are worth emphasizing. 

(i) Interaction matrix element Vk_k' is a-indcpendcnt, 
because 0k (x) — (f>'^\'x.). 

(ii) For a 7^ 0, since e_k 7^ Ck, we have ^_k ^ ^k, 
i?_k 7^ -Ek and /_k 7^ /k- However, we still have U-u = 
C/k, ^'i''k ^ -^k*^ ^^"^ A_k = Ak, as one can see from 



Eqs. UM, UM and ([20]). 

(iii) Fermi energy , relative to which energies such 
as ^k and i?k are measured, is a local quantity, i.e., e^^ = 
ep(a(x')). This can be understood as follows. In the 
superconducting ground state, all electrons are paired. 
For each pair of electrons, we have (ek + e-k) 7^ 2e^"'' for 
a ^ 0. Thus, we expect 7^ e^^"' for a 7^ 0. We further 
expect that > e^°) for a 7^ in the superconducting 
state, where the increase in the Fermi energy is due to a 
flow of the paired electrons. 

For an electronic energy spectrum of the form of Eq. 
(I6|), Fermi energy is determined by = {e^^ + 
e_k )/2, which gives 



2mc2 



(21) 



With the help of Eqs. ^ and ([2T|) . we also have 

Ck + C-k _ ^(0) 
o ~ ^k 



and 



where 



^li^^^k.a, 

2 mc 



AO) _ (0) _ (0) 

Sk ^ ^k ■ 



(22) 



(23) 



(24) 



From the above-described diagonalization of Hamilto- 
nian H, we see the following: Electronic states (or quasi- 
particles) in the superconducting state are each charac- 
terized by a wave- vector k and a spin a, and are in one-to- 
one correspondence with those in the normal state (this 
point is the same as in the case of a = 0)J^ The su- 
perconducting state is different from the normal state 
because (i) there exists an energy gap parameter Ak for 
an electronic excitation in the superconducting state; and 
(ii) Fermi energy is vector potential a dependent in 
the superconducting state. Both the existence of Ak and 
the a-dependence of e^, originate from the pairing of (kf) 
and (— k|) electrons. 
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D. Free energy density F 

From diagonalized Hamiltonian H [Eq. pO)) ]. the fol- 
lowing expression for free energy density F at location x 
(we can drop the prime in x' hereafter) in the supercon- 
ductor can be derived: 



(25) 



F = YXu^-kBT\n{\ + e-''-/''-^]{l + e-^--"'-^' 

k 

ne? 2 (Vxa)^ 



Stt 



The first of the three terms in the above expression 



comes from — /crTIu 



Tr (e-^/'^^^) 



which is the usual 



statistical electronic free energy density.— The second 
term comes from n{ep — t^^^), where n is the density 



of electrons, and we have used Eq. (pij) . This term is 
added to F because electronic energies in the expression 
for Hamiltonian H are measured relative to . The third 
term is the magnetic field energy density. 

Note that the expression for F involves |Ak| (through 
t/k, -Ek and -E-k)- As we will see in the next subsection, 
|Ak| is a function of k, T and a, i.e., |Ak| = |Ak(T,a)|, 



and the function |Ak(T, a)| is determined by the self- 
consistency of the diagonalization of H . Thus, we see 
that F is a function of T and a, i.e., F — F{T, a). [As can 
be shown, free energy density F, excluding the magnetic 
field energy density, becomes a-independent when | Ak| = 
for all electronic states.! 



E. First equation: Implicit solution for |Ak(T, a)| 

The following equation is derived as a condition for the 
self-consistency of the diagonalization of Hamiltonian H: 



1 - /k - f- 



is) 



= independent of T and a . 



(26) 



For a = 0, this equation was previously derived by the 
author in Ref. H. We present the details of the derivation 
of this equation for a 7^ 0, which is similar to that for 
a = 0, in Appendix [Al 

By using 1 — 2/k = tanh(£'k/2fcBr) and the condition 
that |Ak| = at (T, a) = (TcO), we can express the 
above equation as 



lAkP 



tanh 



tanh 



' + |Ak|2-— k-a I /2kBT 
mc 



tanh 



(27) 



This equation is an implicit solution for | Ak | as a func- 
tion of k, T and a for given Tc. 

Note that interaction 14;. k' does not appear in Eq. 
(|27)) . Instead, critical temperature Tc is involved through 
the condition that |Ak| = at (T, a) = (T^O). Namely, 
|Ak| depends on Vk,k' only through T^. 

Critical temperature Tc and phase 0^ of Ak (i.e., 
Ak = lAkle'^"*) are determined by solving the eigenvalue 
problem 



A,< 



k' 



Vk.l 



tanh(|^^?|/2fcBrc) 



(28) 



for given interaction Vk.k' and electronic energy spectrum 
. Equation (^5]) is the linearized form of Eq. (^0]) in 
the limit of (T, a) {Tc, 0). A complete solution for Ak 
is therefore a combination of the solution for | Ak| of Eq. 
(P7)) and the solutions for Tc and 9^ of Eq. (Related 
discussions are also given in Refs. 9 and lift) 



F. Second equation: Relation between j and a 

Having obtained an expression for F{T, a) [i.e., the ex- 
pression given by Eq. (pS)) . with |Ak(T, a)| being imphc- 
itly given by Eq. (|27p ]. we are now ready to consider 
determination of a(x). In thermodynamic equilibrium, 
the overall free energy, given by the volume integral of 
F (r, a(x)), must be stationary with respect to arbitrary 
variation of a(x). This variational problem leads to 



he 



mc m ^ — ' 



(29) 



where j is the electrical current density, and we have used 
the relations = Vxb — VxVxa, with b = V x a 
being the magnetic flux density. 

The first term on the right-hand side of Eq. (^5]) repre- 
sents a uniform flow of all the electrons, which originates 
from the pairing of (kf ) and (— k|) electrons. The second 
term is due to quasi-particle excitations and de-paired 
electrons, and tends to cancel the first term. When all 
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TABLE I: Units for forming dimensionless quantities 



Quantity 



Unit 



Length Ao 

Temperature Tc 

Energy fcsTc 

Magnetic field Hco 

Vector potential XqHco 

Electrical current density (c/47r)_H'co/Ao 

Energy density H^q/Stv 
Effective superconducting electron density n 



the electrons are in the superconducting ground state, the 
second term is zero. On the other hand, when |Ak| — 
for all the electrons (i.e., when the superconductor is in 
the normal state), the two terms cancel each other, and 
we have j = 0. 

Equation (|29p was also derived by BCSJ- However, the 
BCS derivation of Eq. ([^ is based on a linear-response 
approach (i.e., vector potential a is treated as a small per- 
turbation), and thus, is valid only in the low-field limit. 

Note that Eq. ((29|) is not gauge invariant (this feature 
is the same for the London cquation^^). This is because 
the pairing of (k f) and (— k |) electrons fixes the total 
(canonical) momentum of the electrons, and thereby also 
the gauge of the vector potential. Namely, the theory 
became not-gauge- invariant at the point when the pairing 
assumption [Eq. ^] was made. For the same reason, the 



expression for -F(T, a) [Eq. ([25)) ] and the implicit solution 
for |Ak(T, a)| [Eq. (^7)) ] also are not gauge invariant. 
Since, as can be shown, dF/da > (which means that 
a larger value of a is energetically less favorable), we see 
that the gauge of vector potential a in the present theory 
is such that a — > deep inside a bulk superconductor, 
i.e., the London gaugei^i 



G. 



Dimensionless units 



It is convenient to introduce a set of units so that 
physical quantities involved in the theory become di- 
mensionless. The units that we choose to use are listed 
in Table 1, where Aq is the magnetic field penetration 
depth at T = in the limit of zero magnetic field, 
and satisfies Aq — mc? / ATine^ ; and Hco is the thermo- 
dynamic critical magnetic field at T = 0, and satisfies 
i/2^/87r = {-k"^ /&){kBTcfNQ, where iVo is the density 
of states at the Fermi level (A^o = mkj^ /2Tr^h^ for free 
electrons^). The result for Hco was previously obtained 
by the author in Ref. d. 

We also make the substitution (I/Stt'^) Jd^k, 

and, as usual, assume that e^^-* ^ fc^Tc so that the sub- 
stitution de*^") — > and the approximation 
k • a = fcazk — kpttZ]^ hold, where = cosak, with ak 
being the angle between k and a. 

By using the units shown in Table 1 and the above- 
mentioned assumption and approximation, Eq. (j25p be- 
comes 



J 



F = dzk [c/k - Tin (l + e-^-/^) (l + er^--l^)\ + + (Vxa) 



with 



and 



(30) 



and Eqs. (HZD and (HH) become 
2 



(«r)- 



ci.=?r-v'(«f) +iAi<p+iAki 



tanh 



(lefl/2) 



i?±k-^(ef) +|Ak|2±^azk; 



tanh 




azk /2r 



tanh 



+|Ak|^-^azk j/2T 



r 



tanh 



(ieV2) 



(31) 



and 



J = -n^a, 



respectively, where the "effective superconducting elec- 



(32) 
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tron density" 



n„ = 1 



3^2 



na 



/ d^i"^ / dzk^k(/-k-/k) (33) 

Jo JQ 



magnitude of vector potential a, i.e., Ug = ns{T,a) and 
F, = F,{T,a). 

We analyze the functions |Ak(T, a)|, ns{T,a) and 
Fc(T,a) in the next subsection. 



with 



/±k 



Note that Eq. 



+|AkP±^a.k )/T 



(34) 



is the same as the London 



equation, except that our n^, given by Eq. p3p . is both 
T- and a-dependent, whereas ris is only T-dependent in 
the London theory. 

The condensation energy density is defined as 



Fr=FL 



F' 

^ s J 



(35) 



where F!. — i^ — (V x a)^ is the superconducting state free 
energy density excluding the magnetic field energy den- 
sity; and i^/j is the normal state counterpart of F^'. The 
quantity Fc can be used as a measure of the difference 
between the normal and superconducting states. 
As can be shown, the following relation exists: 



2a da 
1 dF, 
2a da 



(36) 
(37) 



where the second expression holds because F^^ is a- 
independent. 

For an isotropic superconductor (as in the present 
case) , Us and Fc are functions of temperature T and the 



H. 



|Ak(r,a)i, n,{T,a), and F,{T, 



We solve Eq. ([3T|) to obtain |Ak(T, a)| by using an 
iterative method. [Note that the variables (k, T, a) for 
the function |Ak(T, a)| appear in Eq. pTjl in the forms 

of d'^^"''!, T, azk); and remember that azk — acosckk is 
the component of a along k .] When solving Eq. ([3T|) , it 
is important to note the following: 

(i) We define temperature T*^^ as such that Eq. (PT|) 

has no |Ak| > solution for given l^^*^^], T and azk if 

T > T*^. Note that is a function of \^^^\ and az^. 
Only for a = is T*^ = Tc the same for all the electronic 
states. 

(ii) Similarly, we define vector potential magnitude a*j^ 
as such that Eq. ((3T|) has no |Ak| > solution for given 
[^^"■'l, T and azk if a > a*j^. Note that a*-^ is a function 

of |^('"|,Tandzk. 

(iii) Depending on |^|^ '|, T and az^, a |Ak| > solu- 



tion of Eq. f3T|) may be an unstable solution, because 
the normal state solution |Ak| may be energetically 
more favorable. 

Point (iii) can be understood as follows. Note that 

for given l^^^''], T and azk, there are always two possible 
solutions for |Ak| if T < or a < a*i^: the |Ak| > 
solution of Eq. (jST]) and the |Ak| = solution. The free 
energy associated with a pair of (k|,— k|) excitations 
with lAkl > is 



-f'sk - ^k - -^k 



tanh(|^f 1/2) 
2|ef I 



1 + e 



l-l-e 



(38) 



where E^'^ = \/ {^^^^^ + l^kP- For |Ak| = 0, it becomes 



(39) 



Both Fsk and F„k are functions of j^j^. |, T and azk- For 

given 1^^"' I' ^ ^'^'^ '^-^k, the | Ak| > solution is the stable 
solution if Fsu < Fn^. Otherwise, the |Ak| = solution 
is the stable solution. 

We define temperature T^k and vector potential mag- 
nitude Ock as such that the |Ak| — solution becomes 
the stable solution for T > Tck or a > Ock- By definition. 



we have T^k < r*k ^'^'^ "^ck < a*k- Note that a pair of 
(kf , — k^) electrons becomes de-paired (having |Ak| = 0) 
for T > Tck or a > Ock- 

Numerical results for Tck and Cck are shown in Fig. 1 
[Figs. l(a)-l(c)]. Figure 1(a) shows T^k versus l^^^''] for 

different values of azk; Fig. 1(b) shows cick-Zk versus |^^°^ | 
for different values of T; and Fig. 1(c) shows Tck versus 




FIG. 1: (a) Tck versus | for different values of azi^; (b) 

cick^k versus |^^°'| for different values of T; (c) T^k versus 

azk for different values of jCk"'!- The dotted curves show the 
corresponding results for T*^^ or a*^^Zk. 



FIG. 2: |Ak| versus \ for different values of T and azk; 
7ra2k/\/2Aoi?cO = 0, 0.8, and 1.2 in (a), (b), and (c), re- 
spectively. The values of T/Tc are indicated on the curves. 
Unstable solutions are shown by dotted curves. 



az]i for diflFerent values of |Ck I- The dotted curves in 
Fig. 1 show corresponding results for T*^^ and a*i^Zk. 

As shown in Figs. 1(a) and 1(b), both Tck and ac^zw 
are increasing functions of |^^*'''|, meaning that pairs of 
electrons with smaller j^^"^ | become de-paired at lower 

values of T and/or az^ than those with larger \^^^\- Fig- 
ure 1(c) shows that Tck is a decreasing function of azk, 
which also means that Ock^k is a decreasing function of 
T. 

For constant l^k*^], the values of Tck and Cck are low- 
est when k is parallel to a, because a appears in Eq. 
pip only in the form of azk, which is largest (there- 



fore, most effective as a de-pairing force) when Zk = 1. 
For two electronic states with wave-vectors k' and k, re- 
spectively, if |^^?^| = \^^^\ and k is parallel to a, then, 
the following relations hold: Tck'(a) = Tckj|a(fl^k'); and 
o-ck'iT) = ack||a(7')/^k'- 

Note that T^k.p\\a.{a) is the lowest value of Tck(a), i.e., 
T'ckFlla(a) = Tck,min(a)- Similarly, acu^^i^iT) is the low- 
est value of ack(T), i.e., acky||a(r) ack,min(r). The 
hmction Tckp||a(a) [or its inverse function aciip\\a.iT)], 

which is shown by the |'Ck'''l = curve in Fig. 1(c) 
(for Zk = 1), is of particular importance. When T < 
Tckp||a(Q.) [or a < ackp||a(^)] in a region in a supercon- 
ductor [remember that a = a(x) is a location-dependent 




FIG. 3: |Ak(r,a2k)j/lAk(0,0)j versus T for different values 
of 14^°' I and azk; l^^^l/ksTc = and 1.0 in (a) and (b), re- 
spectively. The values of Tvaz]i/\/^\oHco are indicated on the 
curves. Unstable solutions are shown by the dotted curves. 



FIG. 4: |Ak(r,a2k)|/|Ak(0,0)i versus azk for different val- 
ues of and T; l^^^l/ksTc = and 1.0 in (a) and (b), 
respectively. The values of T/Tc are indicated on the curves. 
Unstable solutions are shown by the dotted curves. 



quantity], this region is in a "all-paired state," in which 
|Ak| > for all the electrons. On the other hand, 
when T > rcky||a(a) [or a > ackp||a(r)] in a region 
in a superconductor, this region is in a "partly-paired 
state," in which electrons with Tc]i{az]s) in the range 
rckF||a(a) < Tci^iazi^) < T [or with ack(T) in the 
range ac\^p\\ei{T) < ack{T) < a] become de-paired (hav- 
ing |Ak| — 0), while electrons with Tck(azk) > T [or 
ack.{T) > a] remain paired (having |Ak| > 0). 

Numerical results for | Ak| versus |^^°^ | for different val- 
ues of T and azk are shown in Fig. 2 [Figs. 2(a)-2(c)]. 
Figure 2(a) shows the case of a = 0. In this case, |Ak| is 

a monotonic decreasing function of [^k*^^! for all temper- 
atures below T^; and |Ak| vanishes at the same temper- 
ature Tc for all values of [Ck*^!- 

Figure 2(b) shows an example of the case of < a < 
ackj?||a(0) [note that Tracks ||a(0)/V2 = 1 (in dimension- 
less units)]. In this case, the | Ak| -versus- curve for 
T = is the same as in the case of a = 0. However, 
as T increases, |Ak| for a smaller l^^^-*! is more strongly 
suppressed, and decreases faster, so that the |Ak|-versus- 
l^k*^ I curve eventually becomes non-monotonic, with a 
maximum located away from [^^"'l =0. As T increases 
further, |Ak| for a smaller |^k"''l vanishes at a lower tem- 



perature [namely, Tck is smaller for smaller |Ck li a fea- 
ture that is also shown by Fig. 1(a)]. The |Ak|-versus- 
[^^"■'l curve then has two parts: a |Ak| = part for low 
energies, for which rck(a2:k) < T [or ack(T) < a], and a 
|Ak| > part for higher energies, for which Tck(azk) > T 
[or ack(T) > a]. 
Figure 2(c) shows an example of the case of a > 

ackF||a(0). In this case, the |Ak [-versus- [^^"^ | curve has 
a |Ak| = part even for T = 0. Namely, at T = 0, 
|Ak| = for those electronic states with Ock(O) < a. The 
vertical rises in the |Ak [-versus- 1^^^"* | curves for low tem- 
peratures (i.e., the curves for T = 0; 0.3 and 0.5) in Fig. 
2(c) indicate discontinuities. The doted curves in Fig. 
2(c) show corresponding unstable solutions of Eq. ([31]) 
for T in the range Tck < T < T*-^^ (or for a in the range 
flck < a < a^k)- 

Figures 3 and 4 show, respectively, the T-dependence 
and azk-dependence of |Ak|. The dotted curves in Figs 3 
and 4 show corresponding unstable solutions of Eq. ([31]) 
in the range Tck < T < T*^^ or Ock < a < a*^. As shown 
in the figures, |Ak| is a monotonic decreasing function 
of T and azk, except that at T = 0, |Ak| > is a con- 
stant for a < ack(O). Note that, at T = Tck ot a — flck, 
|Ak| may become zero continuously or discontinuously, 
depending on [^k*^^!, T and az^. The vertical drops in 
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FIG. 5: (a) "Effective superconducting electron density" Us 
versus T for different values of a as indicated on the curves; 
(b) Ua versus a for different values of T as indicated on the 
curves. The dotted curve in (a) shows Miihlschlegel's resullii^ 
for Ag/Ai(r), where \l{T) = limif„_o \{T,Ha). 

some of the curves shown in Figs. 3 and 4 indicate dis- 
continuities. 

Based on numerical solutions for |Ak(T, a)|, we can 
compute ns{T,a) and Fc{T,a) straightforwardly by nu- 
merically carrying out the integrals involved in the ex- 
pressions for ns and F^. Numerical results for ns{T,a) 
and Fc{T,a) are shown in Figs. 5 and 6, respectively. 

The dotted curve in Fig. 5(a) shows Miihlschlegel's nu- 
merical result^i^ for Ao/A|(r), which is based on a Hnear- 
response approach,i and therefore, is valid only in the 
low-field limit, i.e., AL(r) = limj^^^o A(T, iJ^). [More 
discussion on magnetic field penetration depth A(T, Ha) 
is given below in Sec. IIIIBI ] Note that Ag/A|(r) cor- 
responds to our ns{T,a)/n for a = 0. The quantita- 
tive difference between Miihlschlegel's Ag/A|(r) and our 
ns(r, 0)/?!, as shown in Fig. 5(a), is due to the use of the 
cut-off approximation^ for solving the energy gap equa- 
tion in Miihlschlegel's work. 

The Fc-versus-T curve for a = in Fig. 6(a) shows 
Fc(r, 0), which is the same as HI{T)/Stt. This result 
was also previously obtained by the author in Ref. 9. 

As shown in Figs. 5 and 6, ns(T, a) and Fc{T^ a) both 
are monotonic decreasing functions of T and a, except 
that Us — 1 is a constant for a < ackj?||a(0) (where 
7rack,.||a(0)/\/2 = 1) at T = [Fig. 5(b)]. Note that. 



FIG. 6: (a) Condensation energy density Fc versus T for dif- 
ferent values of a as indicated on the curves; (b) Fc versus a 
for different values of T as indicated on the curves. 

for T = and a < acky||a(0), we have Fc = 1 — [which 
can be derived either from the expression for Fc, or by 
solving Eq. ^ for = 1 and -Fc(0,0) = 1]. This is 
shown in Fig. 6(b), where we can see that the -Fc-versus-a 
curve for T = is parabolic for < na/ < 1. Namely, 
in this case, all the electrons are paired and no thermal 
excitations exist, so that the decrease in Fc (i.e., the —a^ 
term) is entirely due to the kinetic energy associated with 
the uniform flow of all the electrons. At T = 0, as a in- 
creases further so that a > ackjr||a(0) (i.e., 7ra/-\/2 > 1), 
de-paired electrons (having |Ak| = 0) begin to appear, 
and we have < 1. 

The present analysis provides a possible explanation 
for the experimentally observed non-vanishing Knight 
shifts in the superconducting state near T — (for ex- 
ample, Refs. [ij and ITsI). Namely, at T = 0, electrons 
in a superconducting sample are all paired only when 
the applied magnetic field is weak so that we have a < 
ackjj||a(0) everywhere in the sample; this gives zero spin- 
polarization, and therefore, zero Knight shift (which is 
proportional to the density of spin polarization^^) . How- 
ever, we note that the applied magnetic fields that were 
used for the Knight shift measurement a^^'^^ are compa- 
rable to the thermodynamic critical magnetic fields of 
the samples. Therefore, it is likely that we actually had 

> Q-ckFllalO) over a significant portion of the sample, 
where a finite fraction of the electrons were de-paired and 
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spin-polarized, giving rise to a non-zero Knight shift. The 
previous theoretical prediction of a zero Knight shift in 
the superconducting state for T ^ by Yosida— is valid 
only for weak magnetic fields. Although the present the- 
ory is developed for highly-local superconductors, we ex- 
pect qualitative conclusions of the theory, including the 
prediction for a non-vanishing Knight shift near T = 
for a not-so-weak magnetic field, to be valid also for non- 
local superconductors. 

III. APPLICATIONS 

We present in this section a few examples of applica- 
tion of the theory. We consider the case of a semi-infinite 
superconductor in an applied magnetic field Ha parallel 
to the surface of the superconductor and the case of an 
isolated vortex in an infinite superconductor. We deter- 
mine, in each case, spatial variations of vector potential 
a, magnetic flux density b, electrical current density j, 
energy gap parameter amplitude |Ak|, "effective super- 
conducting electron density" Ug and condensation energy 
density Fc- We also calculate magnetic field penetration 
depth A(T, Ha) and lower critical magnetic field Hci{T). 

A. Semi-infinite superconductor 

We consider a semi-infinite superconductor in an ap- 
plied magnetic field Hq parallel to the surface of the su- 
perconductor. Let the superconductor occupy the half 
space a; > and Hq be applied along the z-axis. In terms 
of the Cartesian coordinates {x, y, z) and the unit vec- 
tors (x,y,z), we can write b — b{x)z, j = i{x)y, and 
a = —a{x)y. Then, Eq. ([5^ and the relation b = Vxa 
become 

h'{x) = -ns(x)a{x) (40) 

and 

a'{x) = ~h{x) , (41) 

respectively, where a "prime" indicates a derivative with 
respect to x; ns{x) — ns{a{x)) is given by Eq. ((33)l . 
where |Ak(a;)| = |Ak(a(x))| is determined by Eq. ([3T|) : 
and the relation j{x) — —b'{x) has been used. 

This is a non-linear second-order boundary- value prob- 
lem with boundary conditions 

b{0) = -a'(0) = Ha (42) 

and 

a(oo) = b{oo) = 0. (43) 

Note that, for the convenience of numerical calcula- 
tion, we have expressed this second-order boundary- value 
problem as a system of two first-order differential equa- 
tions. The numerical method for solving this boundary- 
value problem is explained in Appendix IB] 



Numerical results for a{x), b{x), j{x), ns{x), 
|Akjj||a(2;)| and Fc{x) near the surface of the supercon- 
ductor for different values of Ha are shown in Fig. 7 [Figs. 
7(a)-7(c)] and Fig. 8 [Figs. 8(a)-8(c)] for T = and 0.6, 
respectively. As an example of |Ak(a;)|, |Aky||a(2;)| is 
shown in the figures. Here Akp,||a is the energy gap pa- 
rameter of an electronic state on the Fermi surface with 
wave vector k^? parallel to a, of which the value of Tck(a), 
or ac\t.{T), is the lowest among all the electronic states, 
i-e., Tckj.||a(a) = rck,min(a), and ackF||a(r) = ack,min(T), 
as we mentioned earlier. 

Figures 7(a) and 8(a) show examples of the case where 
Ha is low so that a < a^ky ||a(T') at the surface of the su- 
perconductor. In this case, no de-paired electrons exist 
in the superconductor, i.e., |Ak| > for all the elec- 
trons in the superconductor. At T = [Fig. 7(a)], |Ak| 
and ris are x-independent [because |Ak| and Ug are a- 
independent for a < ackp||a at T = 0, as shown in Figs. 
4 and 5, respectively], and therefore, our results for a{x), 
b{x), j{x) and Fc{x) are the same as those of the London 
theory. Namely, in this case, our Eq. (f32|) is the same 
as the London equation, which gives linear relations be- 
tween a{x), b{x) and j{x), i.e., j{x) cx b(x) cx a{x). For 
T > [Fig. 8(a)], ]Ak| and Ug become cc-dependent near 
the surface of the superconductor, and therefore, our re- 
sults for a{x), b{x), j{x) and Fc{x) near the surface of 
the superconductor are no longer the same as those of 
the London theory. 

Figures 7(b), 7(c), 8(b) and 8(c) show examples of the 
case where Ha has been increased to a point such that 

> o,cU.p\\a.iT) at the surface of the sample. Let xq be 
the coordinate of the location at which a = ackF||a(7")- 
We have )Aky||a(2;)| = in the region < x < xq, 
and |Akp||a(a;)| > in the region x > xo- The region 
< a: < j;o is in a "partly-paired state," where paired 
electrons (having ]Ak| > 0) and de-paired electrons (hav- 
ing ]Ak| = 0) co-exist. The region x > xq remains "all- 
paired," where ]Ak| > for all the electrons. Our results 
for a{x), b(x), j(x) and Fc{x) in the region < x < xq 
are significantly different from those of the London the- 
ory, especially when Ha is high [Figs. 7(c) and 8(c)]. 
In particular, j(x) is no longer a monotonic decreasing 
function of x, having a maximum located near xo- 

B. Magnetic field penetration deptli \{T,Ha) 

The solution for a{x) obtained in Sec. IIII Al for the case 
of a semi-infinite superconductor can be used to calculate 
magnetic field penetration depth X{T,Ha), which, for a 
semi-infinite superconductor, is defined a&i^ 

A = -L rb{^)dx^^, (44) 

ila Jo ila 

where a(0) is the value of a{x) at the surface of the su- 
perconductor. 

Numerical results for A(T, Ha) are shown as Ap/A^ ver- 
sus T for different values of Ha in Fig. 9(a), and Ap/A^ 
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FIG. 7: Variations of a{x), b{x), j{x), na{x), |Akjj||a(a;)|, and 
Fc{x) near the surface at T = and different values of Ha'- 
Ha/H^o = 0.4 (a), 0.6 (b), and 0.8 (c). The quantities a{x), 
b{x), j{x), ns{x), I Akjj||a(a;)|i and Fc{x) are measured in units 
of XoHa, Ha, Ha/Xo, ns{oo), |Akj^(oo)!, and -Fc(oo), respec- 
tively. 



FIG. 8: Variations of a{x), b{x), j{x), na{x), |Akp||a(a;)|, and 
Fc{x) near the surface at T/Tc = 0.6 and different values of 
Ha-. Ha/Hco = 0.3 (a), 0.4 (b), and 0.5 (c) [Ha/HdT) = 0.49 
(a), 0.65 (b), and 0.81 (c)]. The quantities a{x), 6(x), 
j{x), ns{x), |Akjj||a(a;)|, and Fc{x) are measured in units of 
XL{T)Ha, Ha, Ha/XL{T), n^cx)), |Ak^(cxj)|, and F,{oo), re- 
spectively, where Xl{T) — limHa~>o X{T, Ha)- 



versus Ha for different values of T in Fig. 9(b). 

The dotted curve in Fig. 9(a) shows Miihlschlegel's 
numerical result^"^ for A§/A|^(T). As we mentioned 
earlier, the Muhlschlcgel's result is based on a linear- 
response approach, and corresponds to our Us (T, a) / n 
for a = 0, which is the same as X^/\^{T,Ha) for 
Ha = 0, i.e., XllXliT) = lim,^on,(r,a)/n = 
\miHa^o^o/^^{T,Ha). The quantitative difference be- 
tween Miihlschlegel's Ag/A|(r) and our Us (T, a)/n for 
a = 0, or Xl/X^lT, Ha) for Ha ~ 0, is due to the use of 
the cut-off approximation^ in Miihlschlegel's work. 

As shown in Fig. 9(b), for Ha up to about Hc{T)/2, 



Aq/A^ is iJjj-independent or nearly iJa-independent for 
low temperatures, and only weakly iJa-dependent for 
higher temperatures. However, for Ha above about 
Hc{T)/2, Ag/A^ drops rapidly as Ha increases. Note that 
the case of Ha below about Hc{T)/2 corresponds to the 
case shown in Figs. 7(a) and 8(a), where a < acki?||a(2^) 
at the surface of the superconductor so that no de-paired 
electrons exist in the superconductor, whereas the case of 
Ha above about Hc{T)/2 corresponds to the case shown 
in Figs. 7(b), 7(c), 8(b) and 8(c), where a > a^k^\i^{T) 
at the surface of the superconductor so that de-paired 
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FIG. 9: (a) Ag/A^ versus T for different values of Ha as indi- 
cated on the curves; (b) Ag/A^ versus Ha for different values 
of T as indicated on the curves. The dotted curve in (a) shows 
the Miihlschlegel's result^ for Aq/AKt). 



electrons exist in the region near the surface of the su- 
perconductor. 

As shown in Figs. 9(a) and 9(b), Aq/A^ — > when 
Ha Hc{T), or T ^ Tc{Ha) [here TdHa) is the in- 
verse function of Hc{T)]. However, we should note that 
this is not true for a type-II superconductor, for which 
superconductivity is completely suppressed only at the 
upper critical magnetic field Hc2, which is usually much 
higher than He for a highly-local superconductor. This 
is because the definition for A, i.e., Eq. ([^ . is valid 
only when the superconductor is in the Meissner state. 
In the mixed state of a type-II superconductor, in which 
vortices exist, this definition for A is no longer valid. Fur- 
ther, we note that the definition of Eq. (|44p is valid only 
for bulk samples with dimensions much larger than the 
magnetic field penetration depth. 



C. Isolated vortex 

We consider an isolated vortex in an infinite supercon- 
ductor. Let the vortex be centered on the z-axis. In 
terms of the cylindrical coordinates (r, 0, z) and the unit 
vectors (f, 0,z), we can write b — 6(r)z, j = j{r)(f>, and 
a = —a{r)(j). Then, Eq. and the relation b = Vxa 



become 



and 



6'(r) = — ns(r)a(r) 



a'{r) ^ -b{r) 



a{r) 



(45) 



(46) 



respectively, where a "prime" indicates a derivative with 
respect to r; ris(r) ~ 7is(a{r)) is given by Eq. ((33)) . where 
|Ak(r)| = |Ak(a(r))| is determined by Eq. ([31]); and the 
relation j{r) — —b'{r) has been used. 

This is a non-linear second-order boundary- value prob- 
lem with boundary conditions 



b{oo) = a(oo) = 



and 



drrb{r) 



(47) 



(48) 



The last boundary condition comes from fiux quanti- 
zation, i.e., $0 — 27r drrb{r) in conventional units, 
where $o is the flux quantum. The parameter k is de- 
fined as 



2ttXIH,q/^q. 



(49) 



Note the difference between the present definition for 
K and the one in the Ginzburg-landau (GL) theoryi^^ 
^(GL) _ 2-s/27rA^i?c/<i'o- Besides the difference between 
X^Hco for K of Eq. (gH) and X^H^ for there is an 

extra factor in the expression for k'^'^^\ In the GL 
theory^ k^*^^) can also be expressed as a ratio between 
the magnetic field penetration depth and the coherence 
length. A similar expression for k does not exist in the 
present theory, because the coherence effect (or the non- 
local effect) is not accounted for in the present theory. 

For the convenience of numerical calculation, the 
above-described second-order boundary-value problem 
has been expressed as a system of two first-order differ- 
ential equations. The numerical method for solving this 
boundary- value problem is explained in Appendix [Cl 

Numerical results for a (r), 6(r), j(r), ns(r), |Aky||a(r)| 
and Fc{r) near the vortex core are shown in Fig. 10 for 
K = 5 at several different temperatures as indicated in 
the figure. Here, as an example of |Ak(r)|, |Aky||a('^)| is 
shown in the figure. 

Let tq denote the location at which a — ackF||a(r). We 
have |Ak^||a(r)| = for < r < ro, and |Ak^||a(r)| > 
for r > tq. Paired electrons (with |Ak| > 0) and de- 
paired electrons (with |Ak| = 0) co-exist in the region 
< r < To, whereas no de-paired electrons exist (i.e., 
|Ak| > for all the electrons) in the region r > tq. 

As shown in Fig. 10, as r — > 0, we have a —^ oo, 
Us ^ 0, j ^ and Fc 0. It also appears that b{r), 
j{r) and ns{r) have zero slopes at r = 0, while Fc{r) 
has a finite slope at r = and is nearly linear for small 
r. The electrical current density j{r) has a maximum 
located near tq. 
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FIG. 10: Variations of a(r), fe(r), j(r), n,{r), !Akp||a(r)|, 
and Fc{r) near the vortex core at T/Tc = 0.0 (a), 0.6 (b), 
and 0.9 (c), for k = 5. The quantities a{r), b{r), j{r), 
ns{r), I Akjj ||a(?')| , and Fc{r) are measured in units of XoHcO, 
6(0), cHco/4TT\ons{oo), ns(oo), |Akj^(oo)l, and Fc{oo), re- 
spectively. 



D. Lower critical magnetic field Hci{T) 
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FIG. 11: Hci/Hc versus k at different values of T/Tc as 
indicated on the curves. The two curves for T/Tc = 0.6 and 
0.9 are practically indistinguishable from each other. The 
dotted curve shows the Ginzburg-Landau results for Hci/Hc 
in the high-K limit. 



Numerical results for Hd are shown in Fig. 11 as 
Hci/Hc versus k for several different temperatures. 

The dotted curve in Fig. 11 shows the Ginzburg- 
Landau result for Hd/Hc for high-K superconductorsi*^ 



H. 



(GL) 



1 



(lnK+ 0.85), 

2k 



(hiK^'^-^^ + O.So) 



(51) 



where we have used k^'-^^' — \/2k by replacing X^Hc with 
XqHco in the expression for k'^'^^K 

As shown in Fig. 11, our result for Hd/Hc is only 
weakly T-dependent for low temperatures, and nearly 
T-independent for intermediate and high temperatures, 
as indicated by the fact that the two Hd/ Hc-veisns-K 
curves for T = 0.6 and 0.9 are practically indistinguish- 
able from each other. This feature is to be compared 

(GL) 

with that H^^ /He is T-independent in the Ginzburg- 
Landau theory. Figure 11 also shows that our result for 
Hd/Hc and the Ginzburg-Landau result are quantita- 
tively not very different for k 3> 1. 

For lower values of k, the present theory underesti- 
mates the value of Hd, because the coherence (or non- 
local) effect, which increases the energy associated with 
a vortex and thus leads to a larger Hd , is not accounted 
for in the present theory. 



The numerical solutions obtained in the last subsec- 
tion for an isolated vortex can be used to calculate lower 
critical magnetic field Hd of the superconductor. By 
definition, at Ha — Hd, the Gibbs free energy must 
be the same whether the first vortex exists or not, i.e., 

^(no vortex) ^ ^ (one vortex) ^ ^^.^ COudition 

leads to 

Hd = ^ drr {H^ -Fc + b") . (50) 



IV. SUMMARY 

We have presented a microscopic theory for supercon- 
ductivity in a magnetic field based on a local approxima- 
tion approach. The theory allows microscopic description 
of the suppression of superconductivity by an externally 
applied magnetic field. 

In Sec. ini we presented the details of the theory. The 
main results derived in Sec. [Til include an expression for 
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free energy density F as a function of temperature T 
and vector potential a, and two basic equations of the 
theory: the first is an imphcit solution for energy gap 
parameter amplitude |Ak| as a function of wave vector 
k, temperature T and vector potential a; and the sec- 
ond is a London-like relation between electrical current 
density j and vector potential a, with an "effective su- 
perconducting electron density" that is both T- and 
a-dependent. The two equations allow determination of 
the spatial variations of a and |Ak| in a superconduc- 
tor for given temperature T, applied magnetic field 
and sample geometry. In the low- field limit, the the- 
ory reduces to the local-limit result of BCS. We also nu- 
merically analyzed the functions |Ak(T, a)|, ns{T, a) and 
Fc{T,a) (where Fc is the condensation energy density). 

In Sec. mil as examples, we applied the theory to 
the case of a semi-infinite superconductor in an applied 
magnetic field Ha parallel to the surface of the supercon- 
ductor and the case of an isolated magnetic vortex in an 
infinite superconductor, and determined, for each case, 
spatial variations of quantities such as a, |Ak|, and 
Fc- We also calculated magnetic field penetration depth 
X{T,Ha) and lower critical magnetic field Hci{T). 

An important conclusion of the theory is that, depend- 
ing on temperature T, applied magnetic field and 
sample geometry, a "partly-paired state" can exist in 
which paired electrons (having |Ak| > 0) and de-paired 
electrons (having |Ak| = 0) co-exist. Such a "partly- 
paired state" exists even at T = when Ha is above 
a threshold for a given sample, giving rise to a non- 
vanishing Knight shift in the superconducting state at 
T = for Ha above the threshold. 

Since the non-local effect (or coherence effect) in the 
superconducting state is not accounted for in the present 
theory, we expect the theory to be valid only for highly- 
local superconductors (for which magnetic field penetra- 
tion depth A is much larger than coherence length ^). 
However, when a more complete theory is developed that 
is able to account for the non-local effect (or coherence 
effect), we expect it to reduce to the present theory in 
the local limit. 



APPENDIX A: DERIVATION OF EQ. ([26 



tion, Eq. ([20)) . can be rewritten as 



For a = 0, Eki. ([^5]) was previously derived by the 
author in Ref. |8| (related discussions are also given in 
Refs. i and [13). The derivation of Eq. ([261) for a ^ is 
similar to that for a = 0. We now present the details of 
the derivation of Eq. (^5)1 . 

For convenience, we define 



1 - /k - /- 



2E, 



is) 



(Al) 



Ak - -Xl^kk'Ck'Ak' 



(A2) 



In the presence of an applied magnetic field, we ex- 
pect Ak to be a function of temperature T and vector 
potential a = (ai, 02,03), where ai {i = 1,2,3) are the 
components of a. Let X denote any one of T, ai, 02 
and 03. We operate d/dX on both sides of Eq. (|A2[) to 
obtain 



5Ak 

dX 



. f ^A. 
^ \dX ^ 



dX 



(A3) 



We next multiply both sides of the above equation by 
CkAJ^, and then take summation over k, i.e.. 



aAk 

dX 




f^A 

I dx 



(A4) 



The quantity inside the first pair of parentheses on 
the right-hand side of the above equation equals to AJ^, 
[according to Eq. (|A2[) ] so that the second of the two 
terms on the right-hand side is the same as the term on 
the left-hand side. Thus, we have 



.5Ck 
dX 



= 



We want a |Ak| > solution. Clearly, 

Ck — independent of T and a. 



which is Eq. (|26|) and satisfies 
5Ck 



dX 



0, 



(A5) 



(A6) 



(A7) 



is a solution of Eq. (jA5[) . 

However, since Eq. (|A6|) is not the only possible solu- 
tion of Eq. (|A5P [as one can see, Eq. (|A5|) actually can 
have an infinite number of solutions], we need to justify 
that Eq. (|A6p is the only physical solution. 

Since the diagonalized Hamiltonian, Eq. PH]) . de- 
scribes a set of independent quasi-particle excitations, 
there should be no coupling (except pair correlation) be- 
tween the quasi-particle excitations. Therefore, we ex- 
pect the thermal energy and entropy associated with each 
pair of (kt, — k|) excitations to be 



£k — C^k + /k£^k + /-k^^-k 



(A8) 



and 



^k 



which is a real number. Then, the self-consistency equa- 



-ks [/kln/k + (f-/k)ln(l-/k) 

- /-k In /-k + (1 - /-k) ln(l - /-k)] , (A9) 
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respectively. Similarly, we expect the contribution to the 
electrical current density from each pair of (k f,— k J,) 
excitations to be 



Jk = e (/_k - /k) Vk 



(AlO) 



[here Vk — Kk./m for free electrons, and X^k^^k corre- 
sponds to the second term on the right-hand side of Eq. 



However, as compared to the above expressions for 
£kj ^k, and jj^, those derived from the diagonalized 
Hamiltonian contain additional terms involving dUi^/dX, 
dEk/dX and dfk/dX (where X = T in the case of ek 
or Sk', and X = ai in the case of Letting the sum of 
the additional terms to be zero, one gets Eq. (jA7p . and 
therefore Eq. (|A6p . 



APPENDIX B: NUMERICAL METHOD FOR 
SOLVING THE BOUNDA RY- VA LUE PROBLEM 
OF SEC. IIII Al 



The boundary-value problem of Sec. IIII A| as specified 
by the system of Eqs. (|40|) and (|4T|) . with boundary 
conditions given by Eqs. (|^^ and can be solved by 
using the Runge-Kutta method. 

In order to use the Runge-Kutta method, we first need 
to know values of a and b at one point on the x-axis. 
We know that deep inside the superconductor both a 
and b become small, and Us becomes an x-independent 
constant, so that the London solutions for a and b hold, 
which give the following relation between a and b: 



(Bl) 



Usd. 



Let xq denote the coordinate of such a point located deep 
inside the superconductor. Since both Eqs. ([^0]) and (PT|) 
do not involve coordinate x explicitly, we can assign an 
arbitrary value to xq. We then assign a sufficiently small 
value for a{xo), and obtain b{xo) from the above relation, 
where ns{xo) is obtained from Eq. (|33p for a — (since 
the a-dependence of ris is negligible when a is small) and 
given T. 

Once we know a{xo) and b{xo) at x = xq, we use the 
Runge-Kutta method to compute a{xn) and b{xn) for 
Xn — Xq — nh, where /i is a small positive interval and 
n = 1, 2, • • • , iV, until b(xN) > Ha at xn = — Nh. 

Usually, b{xN) is greater than Ha- If the difference 
between b{x^) and Ha is small, we can simply use x^ as 
coordinate Xg of the surface of the superconductor. Or, 
we can obtain Xs by making a linear interpolation: 



XN-1 



^ Hg -b{xN-i) 
b{xN) - b{xN-i) 



(B2) 



APPENDIX C: NUMERICAL METHOD FOR 
SOLVING THE BOUNDA RY- VA LUE PROBLEM 
OF SEC. IIII Cl 



The boundary-value problem of Sec. IIII Ci as specified 
by the system of Eqs. (|45l) and (|46|) . with boundary 
conditions given by Eqs. P5)l and ([35]), can be solved by 
using the Runge-Kutta method. 

In order to use the Runge-Kutta method, we first need 
to know values of a and b at one point on the r-axis. We 
know that far away from the vortex core, both a and b 
become small and Ug becomes a constant. For a constant 
Us, Eqs. (Uni and can be solved analytically (see, 
for example, Ref. l21))): 



C 



and 



b = CKo{r^), 



(Cl) 



(C2) 



where C is a constant to be determined, and Kn(x) are 
the modified Bessel functions of the second kind. 

Let To denote the r-coordinate of such a point located 
far away from r = 0, the center of the vortex. We choose 
a value for vq that is sufficiently large, and guess an ini- 
tial value for constant C, say Cq, and obtain a{ro) and 
b{ro) from Eqs. (|C1[) and (jC2[) . respectively, where Ug is 
obtained from Eq. p3p for a = (since the a-dependence 
of rig is negligible when a is small) and given T. Once we 
know a(ro) and 6(ro) at r = tq, the Runge-Kutta method 
allows us to compute a(r) and b{r) for any r. 

We then compute total magnetic flux ^E" associated the 
vortex by numerically carrying out the integral on the 
right-hand side of Eq. (pS)) . If total magnetic flux ^ = 
^'o for C = Cq is, for example, greater than flux quantum 
7 we reassign a smaller value for C, say Ci, and repeat 
the computation of a(r), b{r), and 5". 

For the i-th repetition {i > 2), we can assign a value 
for C by making a linear interpolation or extrapolation, 
i.e., 



Ci — Ci-2 + 



^0 - ^»-2 



{Ci-i — Ci-2) 



The computation of a(r), b{r) and 4" is repeated until 
the difference between total magnetic flux VE" and flux 
quantum $0 is within a predetermined range. In practice, 
it usually involves only a few repetitions. 



1 J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev. ^ L. P. Gorkov, Zh. Eksp. Teor. Fiz. 36, 1918 (1959); 37, 
108, 1175 (1957). 1407 (1959) [Sov. Phys. JETP 9, 1364 (1959); 10, 998 



16 



(I960)]. 

^ N. R. Werthamcr, in Superconductivity, edited by R. D. 
Parks (Marcel Dckkcr, New York, 1969). 

* N. N. Bogoliubov, Nuovo Cimento 7, 794 (1958); Zh. 
Eksp.Teor. Fiz. 34, 58 (1958) [Sov. Phys. JETP 7, 41 
(1958)]; J. G. Valatin, Nuovo Cimento 7, 843 (1958). 

^ N. W. Ashcroft and N. D. Mermin, Solid State Physics, 
Holt, Rinchart and Winston, New York (1976). 
A. B. Pippard, Proc. Roy. Soc. (London) A216, 547 
(1953). 

R. P. Feynman, Statistical Mechanics, Benjamin, Reading, 
MA (1972). 

* Z. Hao, Mod. Phys. Lett. B7, 1439 (1993). 

^ Z. Hao, J. Phys. Chem. Solids 57, 1215 (1996). 

Z. Hao, "New interpretation for energy gap A of the cut- 
off approximation in the BCS theory of superconductivity" 
(preprint). 

F. London and H. London, Proc. Roy. Soc. (London) 
A149, 71 (1935); F. London, Superfluids, vol. 1, Dover 

Publications, New York (1961). 

S. D. Conte and C. de Boor, Elementary Numerical Anal- 
ysis: An Algorithmic Approach, 3rd edition, McGraw-Hill, 



New York (1980). 

B. Miihlschlegel, Z. Phys. 155, 313 (1959) [cnglish trans- 
lation in The Theory of Superconductivity, edited by N. N. 
Bogoliubov, Gordon and Breach, Science Publishers, New 
York (1968)]. 

F. Reif, Phys. Rev. 102, 1417 (1956); 106, 208 (1957). 

G. M. Androes and W. D. Knight, Phys. Rev. Lett. 2, 386 
(1959); Phys. Rev. 121, 779 (1961). 

C. P. Slichter, Principles of Magnetic Resonance, Harper 
and Row, New York (1963). 

K. Yosida, Phys. Rev. 110, 769 (1958). 

V. L. Ginzburg and L. D. Landau, J. Exptl. Theoret. Phys. 

(USSR) 20, 1064 (1950) [Emglish translation in Men of 
Physics: L. D. Landau, vol. 1, edited by D. ter Haar, Perg- 
amon Press, New York (1965)]. 

A. L. Fetter and P. C. Hohenberg, in Superconductivity, 
edited by R. D. Parks (Marcel Dekker, New York, 1969). 
C.-R. Hu, Phy. Rev. B 6, 1756 (1972). 
M. Tinkham, Introduction to Superconductivity, McGraw- 
Hill, New York (1975) 



